Time-of-flight mass measurements for nuclear processes in neutron star crusts 
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The location of electron capture heat sources in the crust of accreting neutron stars depends on 
the masses of extremely neutron-rich nuclei. We present first results from a new implementation of 
the time-of-flight technique to measure nuclear masses of rare isotopes at the National Supercon- 
ducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium - nickel range 
were determined simultaneously, improving the accuracy compared to previous data in 12 cases. 
The masses of 61 V, 63 Cr, 66 Mn, and 74 Ni were measured for the first time with mass excesses of 
-30.510(890) MeV, -35.280(650) MeV, -36.900(790) MeV, and -49.210(990) MeV, respectively. 
With the measurement of the 66 Mn mass, the locations of the two dominant electron capture heat 
sources in the outer crust of accreting neutron stars that exhibit superbursts are now experimen- 
tally constrained. We find that the location of the 66 Fe— > 66 Mn electron capture transition occurs 
significantly closer to the surface than previously assumed because our new experimental Q-value 
is 2.1 MeV (2.6a) smaller than predicted by the FRDM mass model. 
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Neutron stars that accrete matter from an orbiting 
low-mass companion star are observed as galactic X-ray 
binaries [ij. A fluid element accreted onto the neutron 
star surface is buried by the continuous accretion of more 
matter, and undergoes a sequence of compositional trans- 
formations driven by nuclear reactions under rising pres- 
sure. Near the surface, at typical depths of a few meters, 
thermonuclear explosions, which are observed as X-ray 
bursts, burn hydrogen and helium into heavier elements 
in the nickel - cadmium range Q. In somewhat deeper 
layers explosive carbon burning is thought to produce 
the occasionally observed superbursts Q , converting the 
ashes of the regular bursts into nuclei in the iron - nickel 
range. Still deeper in the neutron star crust the mat- 
ter undergoes a sequence of electron captures, accom- 
panied at even greater depth by neutron emissions, and 
pycnonuclear fusion reactions [4j-|6|. These nuclear pro- 
cesses, which involve extremely neutron- rich nuclei, heat 
the crust creating a characteristic temperature profile. 

In this letter we report results from an experiment 
at the National Superconducting Cyclotron Laboratory 
(NSCL), where we have produced and measured the 
masses of neutron-rich nuclei using a new implementa- 
tion of the time-of-flight (TOF) technique. We measured 
the mass of sixteen neutron-rich isotopes in the region 
around N = 40, four of which were measured for the 
first time. The results allow us to locate the dominant 
electron capture heat sources in the crust of accreting 



neutron stars and better constrain their strength. The 
masses also provide new information on the onset of de- 
formation near the N = 40 region in neutron-rich nuclei. 

An understanding of crustal heating nuclear processes 
in accreting neutron stars is needed to interpret a number 
of observables. Neutron stars in transiently accreting X- 
ray binaries offer the unique opportunity to directly ob- 
serve the temperature profile of the crust. Some systems 
accrete for many years, sufficiently long to reach ther- 
mal equilibrium. Then accretion stops for many years, 
enabling the observation of the thermally relaxing crust 
over time (for example Ref. Q). The time dependence 
of the cooling curve contains information on crust prop- 
erties such as composition, thermal conductivity, heat 
capacity, neutron superfluidity, and the efficiency of neu- 
trino cooling [H, Q . Interpretation of these observations 
requires reliable nuclear physics to predict the location 
and strength of the nuclear heat sources during the accre- 
tion phase. Crustal heating predictions are also needed 
to understand the recurrence time of superbursts Q , and 
the generation of gravitational waves due to the deforma- 
tions induced by electron capture reactions in the crust 
of the rapidly spinning neutron star, which might be ob- 
servable with future gravitational wave detectors [l(| • 

Which electron captures occur in the outer crust of an 
accreting neutron star depends on the composition syn- 
thesized by thermonuclear burning processes at the neu- 
tron star surface. Model calculations show that for sys- 
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FIG. 1. Integral of the heat deposited in the neutron star 
crust by electron captures as a function of depth (indicated 
through p e ). Shown are results for masses from FRDM 12] 
only (red, dot-dashed), for implementing previously published 
experimental masses (green, solid), and for implementing in 
addition our new mass results (dashed, black). 
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FIG. 2. Black diamonds show fit residuals for the reference 
isotopes as a function of their mass-to-charge ratio ( 51 Sc, 52 Sc, 
54 Ti, " - - - - 

63 Mn, 

75 Zn, 76 Zn, and 79 Zn). Blue circles are the isotopes for which 
we present an improved mass value (see Table [I]). 
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terns exhibiting superbursts, resulting ashes are mainly 
nuclei in the A = 54 - 66 mass range [ll|. These ashes 
serve as the initial composition for the electron capture 
processes. Crust model calculations [1] show that for such 
composition, heat release in the outer crust is dominated 
by two transitions: two-step electron captures on 66 Ni — > 
66 Co 66 F(3j and subsequently on 66 Fe -)• 66 Mn -> 66 Cr 
(corresponding to steep increases in depth-integrated 
heat in Fig.[T]). The location of these transitions is set by 
their electron capture threshold Qec- Because the elec- 
tron chemical potential p c slowly rises with depth and 
the temperature is rather low (kT < 40 keV<§; p e ) the 
transition occurs essentially at a depth where jj, e « Qec- 
Qec = A(Z, N) - A(Z -1,N+1)-E X depends on the 
mass excess A of parent and daughter nuclei, and the ex- 
citation energy of the lowest lying state E x into which the 
capture can occur. Because of nuclear pairing there is a 
strong odd-even staggering of Qec and the threshold for 
the two-step transition is effectively set by the first step, 
the electron capture on 66 Ni and 66 Fe. In both cases E x is 
predicted to be negligible (of the order of 0.1 MeV), and 
therefore it is the nuclear masses of 66 Ni, 66 Co, 66 Fe, and 
66 Mn that determine where the heat is deposited. With 
our first mass measurement of 66 Mn, all these masses are 
now known experimentally. 

TOF mass measurements have been successfully ap- 
plied with different technical approaches at several facil- 
ities for the study of short-lived isotopes [l3|, 14 1. We 



report here a first implementation of the TOF technique 
at the NSCL. Details of the setup and the analysis will 
be presented in an upcoming publication. Neutron-rich 
isotopes were produced by fragmentation of a 86 Kr pri- 
mary beam at 100 MeV/u in a Be target. The frag- 



ments were collected by the A1900 fragment separator 
and transmitted to the S800 spectrometer [l5| through 
a bcamlinc. Two production targets, with thicknesses of 
51 mg/cm 2 and 94 mg/cm 2 , were alternated, keeping the 
magnetic rigidity of the A1900, the beamline, and the 
S800 unchanged, to increase the transmission of parti- 
cles with different mass-to-charge ratios. Thus we obtain 
a sufficient number of calibration nuclei with well-known 
masses. The TOF was measured with newly developed 
fast timing scintillators located at the focal planes of the 
A1900 and S800 resulting in a flight path of 58.7 m. The 
momentum acceptance of the system was 0.5 %, requiring 
a precise relative magnetic rigidity {Bp) measurement of 
each beam particle. This was accomplished with a posi- 
tion sensitive micro-channel plate detector located at a 
dispersive focus of the S800. Detectors at the S800 fo- 
cal plane provided energy loss measurements for particle 
identification and beam tracking information. 

From the simultaneous measurement of magnetic rigid- 
ity, TOF, and atomic charge number (from energy loss) 
for each fully-stripped beam ion, the mass can be deter- 
mined. The measured TOF of each isotope was corrected 
for its dependence on the measured Bp using an empir- 
ical relationship. The resulting relative mass resolution 
was 1.8 x 10~ 4 for the typical case. 

The relation between TOF and m/q of each ion was 
obtained by fitting a 6 parameter calibration function of 
second order in TOF and third order in Z to the mea- 
sured TOFs of 26 reference isotopes of known mass 0- 
Hil ]. The reference masses included two isotopes with 
known low-lying isomers ( 64 Mn and 75 Cu [20]). We have 
confirmed that the unknown population of these isomers 
does not affect the final results, by performing different 
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TABLE I. Mass excess results, in keV, from the present ex- 
periment and from the literature [161 ] (all corresponding to 
measurements at TOFI [2lj]). The fourth column shows the 
counts of each isotope used for our measurement. 
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up to 440 keV 



fits with variations of the 64 Mn and 75 Cu masses that 
account for a range of possible isomeric populations, fol- 
lowing the method in Ref. [l6| • The resulting fit residuals 
show no apparent systematic trends (Fig. f2|) but the \ 2 
per degree of freedom of the fit is larger than one, indi- 
cating the presence of additional systematic errors. To 
estimate the magnitude of these errors we find the addi- 
tional error that normalizes \ 2 P er degree of freedom to 
one when added in quadrature to each m/q calibration 
data point. We find a systematic error of 5.3 keV/q, or 
130 keV for manganese. To determine the error of a mea- 
sured mass, we add in quadrature this systematic error, 
the statistical error, and the calibration error from the 
uncertainties of the fit parameters explicitly calculated 
from the fit covariance matrix. Results are shown in Ta- 
ble |H For the new masses in this work the statistical 
error dominates, with the calibration error contributing 
a significant fraction. 

Figure [3] shows the systematics of two-neutron sepa- 
ration energies (S<2 n ) extended towards more neutron- 
rich nuclei by our measurements. As the fp-shell gets 
filled with neutrons, configurations involving the gg/2 
shell start driving increased deformation near N = 40 for 
Z < 28 nuclei leading to increased binding energies, and 
a change in the slope of the observed two-neutron sepa- 
ration energies as functions of neutron number [23j |. For 
iron, 7-spectroscopy studies have demonstrated a marked 
decrease in the energy of the first 2 + state in even-even 
isotopes starting at N — 38 indicating the onset of in- 




FIG. 3. Two neutron separation energies as a function of neu- 
tron number. S2n values calculated with our masses results 
are shown as filled symbols (red). Open symbols show 
values from the literature 



creased deformation, coinciding with the slope change in 
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24 |25| . A similar effect is observed for chromium 
starting at N = 36, again confirmed by 7-spectroscopy 
studies 2(| 27| . For manganese and iron our new masses 
confirm the continuation of this trend beyond N = 40. 
Interestingly, for vanadium isotopes our new mass of 61 V 
shows the onset of the same effect beyond N = 36. While 
our uncertainty for 61 V is large (890 keV) the deviation 
from the linear trend of the lighter vanadium isotopes is 
about 3.5 MeV, or 4er. This is in line with a comparison of 
experimental /3-decay half-lives with shell model calcula- 
tions, which indicated that for 61 V a pure fp-shell model 
calculation is not adequate, hinting at the onset of the in- 
fluence of the g 9 /2 shell [28|. For the cobalt isotopes, our 
mass results are systematically lower thanprevious mea- 
surements with the TOFI spectrometer 2l|. The same 
trend was observed in recent Penning trap measurements 
for 66 Co and 67 Co 0. 

Using the new masses we carried out model calcula- 
tions of the heating in the crust of an accreting neu- 
tron star, using the same steady state crust model de- 
scribed in detail in The initial composition is the 
ashes of a carbon superburst occurring at a depth where 
fi c w 4 MeV. Shown in Fig.[T]are results using the FRDM 
mass model fl2| . and results when including available 
experimental masses. The FRDM was employed here, 
as well as in previous crust model calculations for 
consistency with the QRPA approach used to calculate 
electron capture rates. The heating is dominated by two 
sources: electron capture on 66 Ni (at fx c « 9.4 MeV with 
FRDM), and on 66 Fe (at fi e « 15.3 MeV with FRDM). 
Previous mass measurements allow one to pinpoint the 
electron capture on 66 Ni, shifting it slightly deeper to 
fi e « 9.6 MeV. The energy this transition releases in the 
crust is smaller because of a smaller odd-even staggering 
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in the experimental Q-values (for 66 Co the experimen- 
tal Q EC = -6.34 ± 0.02 MeV, while the FRDM value is 
Qec — —5.40 MeV). Our new result of the 66 Mn mass 
allows us now to also place the 66 Fe transition based on 
an experimental Qec- Because our new Qec for 66 Fe 
of -13.2 ± 0.8 MeV is 2.1 MeV {2.6a) smaller than the 
FRDM prediction, the 66 Fe transition turns out to occur 
at a much shallower depth, around fi c « 13.2 MeV. 

A possible reason for the discrepancy of the FRDM 
prediction of the Qec for 66 Fe is that the spherical Z=28 
shell gap is too strong in the model preventing the on- 
set of deformation at iV=40. We can also compare our 
new 66 Fe Qec result with Q-values obtained with the 
microscopic mass models HFB-14 and HFB-17 j^. 
HFB-14 predicts Q EC = -12.2 MeV, more than la dif- 
ferent from the experimental value. On the other hand, 
the recent HFB-17 predicts Qec = —13.2 MeV in agree- 
ment with experiment, though the individual masses are 
each discrepant by about 0.9 MeV. However, the error es- 
timates of the Qec theoretical predictions (corrected for 
experimental errors following the procedure in 12j, Eq. 
(6)) in the entire region of interest (even A chains with 
39 < A < 71 and Qec < -3 MeV) for FRDM, HFB-14, 
and HFB-17 are 0.76 MeV, 0.89 MeV, and 0.85 MeV, 
respectively. 

In Fig.Q]our new masses lead to a reduction of the total 
heat produced compared to [1| because the 66 Fe transi- 
tion occurs at a lower fj, e , leading to a smaller fi c — Q 
for the capture on 66 Mn that directly follows. The Qec 
for 66 Mn is still taken from FRDM, potentially leading 
to an unrealistic odd-even staggering. A mass measure- 
ment of 66 Cr is needed to address this issue. In addition, 
the amount of crustal heating depends on the prediction 
of the correct excitation energies for the final states for 
electron captures on 66 Co and 66 Mn, which affects the 
fraction of energy lost by neutrino emission. 

Crustal heating in accreting neutron stars can be 
strongly affected by nuclear structure effects associated 
with sub-shell closures [J]. Our results for A=40 show 
that at the same time, mass models can be particularly 
uncertain in these regions. A similar effect might be ex- 
pected for neutron rich nuclei with A ~100 near A^=60 
[3l| . which would be relevant for models of heating in 
neutron star crusts with a subtantial amount of heavier 
X-ray burst ashes. 

In summary, we have presented the first results from a 
new implementation of the TOF technique at the NSCL 
to measure masses of very neutron-rich nuclei. System- 
atic errors of approximately 130 keV have been achieved, 
and the resolution is sufficient for measurements with a 
few 100 keV accuracy for nuclei where a few 1000 events 
can be produced. Our measurement of the 61 V mass 
indicates increased deformation and the beginning of in- 
fluence of the g 9 / 2 orbital in this isotopic chain. Our first 
measurement of the 66 Mn mass allows us to determine 
the depth of all the major heat sources in the outer crust 



of accreting neutron stars. We find that one of them, 
the electron capture on 66 Fe, occurs at a much shallower 
depth than predicted. The 800 keV uncertainty of the 
new electron capture Q-value for 66 Fe is much smaller 
than typical variations in mass model predictions. Fur- 
thermore, the uncertainty is now experimentally deter- 
mined and can be considered in astrophysical models in 
a quantitative way. With the depth of the transitions 
fixed, the remaining smaller nuclear physics uncertain- 
ties affecting the amount of heat deposited can now be 
taken into account when using models of crustal heating 
to interpret observations. A mass measurement of 66 Cr 
would be helpful to reduce these uncertainties further. 

This work was partially supported by NSF grants 08- 
22648 and PHY 06-06007, and by DFG under contract 
number GE 2183/1-1. 
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